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The recent results in y/s = 2.76 TeV Pb+Pb collisions at the Large Hadron Collider (LHC) re- 
ported by the ALICE collaboration shows that the power-law energy-dependence of charged hadron 
' multiplicity in Pb+Pb collisions is significantly different from p+p collisions. We show that this 

different energy-dependence can be explained by inclusion of a strong angular-ordering in the gluon- 
decay cascade within the Color-Glass-Condensate (or gluon saturation) approach. This effect is 
more important in nucleus-nucleus collisions where the saturation scale is larger than 1 GeV. Our 
prescription gives a good description of the LHC data both in p+p and Pb+Pb collisions. 

I. INTRODUCTION 

. The recent LHC data on hadron production in proton-proton (pp) and nucleus-nucleus (AA) scattering [H-HI shows 
^p H ' that gluon saturation that follows both from the BFKL Pomeron calculus @ and from the Color Glass Condensate 
O r (CGC) approach @-[l(|, gives an adequate description of the high energy scattering in QCD. The model based on 
the gluon saturation was able to predict the hadron production at y/s = 7 TeV [ll[ (see also Ref. [l2j|) and the 
experimental data both for pp [l[ and AA collisions [1, Q confirmed the basic qualitative predictions of this approach 
[Hi [HI- However, the recently reported data from the ALICE collaboration 0] on hadron production in AA collisions 
also demonstrated that we are far away from the high precision quantitative description. For example, the model 
that predicted 7 TeV data for pp scatterings and which also describes HERA and RHIC data, failed to describe the 
multiplicity in AA collisions [15| with the same accuracy. This fact cannot be considered as discouraging since AA 
collisions are more complicated QCD problem and moreover the other model calculations, based on the same ideas, 
were somehow able to describe the data [HI, (HUH)]. Nevertheless, this gives rise to a question that whether despite of 
considerable progress in theory during the past two decades we are ready to give a reliable prediction in the framework 
£N ■ of the high-density QCD. 

In practice, our theoretical description of hadron-hadron and nucleus-nucleus scatterings is based on two main 
I . ingredients: Balitsky-Kovchegov (BK) non-linear equation [10j and k t factorization |l7|,|l8|- However, the BK equation 
• • ' is not complete since it does not take into account the correct (non-perturbative) behaviour at large impact parameters 
and, because of this, it leads to the violation of the Froissart theorem, see Ref. [l^| where this problem discussed in 
details. A practical consequence of this is the fact that we may not be able to guarantee the accuracy better than 
±20% (if not worse) (2pl - [22j | from the application of the BK equation. Having said that, the recent application of the 
BK equation to the description of HERA data looks promising indeed [23[ . On the other hand, the k t factorization 
is not reliable for dense-dense system scatterings [24-26] and, strictly speaking, we cannot apply this factorization 
neither to proton-proton nor to nucleus-nucleus scatterings at midrapidity. Therefore, we are doomed to build models 
trying to get a feedback from the experimental data for a theoretical breakthrough. 

As a first attempt toward understanding of the new LHC data on nucleus-nucleus collisions, one may compare 
the experimental data with the principal qualitative predictions of the gluon saturation. The cornerstone of such 
predictions is the fact that the multiplicity in pp and AA collisions are proportional to Q 2 S cx s x l 2 @, [TH, [13, El, H3> 
where Q s is the saturation scale, s is the center-of-mass energy squared per nucleon pair and A is free parameter 
to be fixed with other experiments like DIS at HERA. This indicates that the energy dependence of multiplicity in 
both pp and AA collisions should be the same, assuming that the atomic number or A dependence of the saturation 
scale is factorizable from energy. This simple property is in accordance with RHIC data [HI, H3, HH . However, the 
new ALICE data shows that multiplicities in pp and AA collisions have a different energy power-law behaviour (see 
Fig. []}. Thus, at first sight, it looks as if that one of the principal feature of high-density QCD is violated. In this 
paper, we will argue that indeed the recent LHC data on AA collisions at yfs — 2.76 TeV has already opened up 
a new QCD regime which requires further theoretical understanding than previously thought. Here, we shall give a 
simple explanation of the different energy behaviour of the hadron multiplicity in AA and pp data at high energy 
based on the gluon saturation (or the CGC) scenario. 

In the CGC approach, the hadron production goes in two stages: production of gluons and subsequently the decay 
of gluon-jet (or mini-jet) into hadrons. Therefore, the multiplicity of the produced hadrons at pseudo-rapidity 77 can 
be calculated as a convolution of these two stages, 
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FIG. 1: The energy behaviour of charged particle pseudo-rapidity per participant pair for central AA and non-singlet diffractive 
pp collisions. The energy dependence can be described based on the saturation picture by s ' 1 for pp and s°' 14j for AA collisions. 
The saturation (CGC) curve for pp collisions is taken from Ref. The saturation curve for the A A collisions was calculated 
from Eq. (O having incorporated the effects of gluon-jet angular ordering which is important when the saturation scale Q s > 1 
GeV, see the text for the details. The total theoretical uncertainties in the saturation model calculation is about 7% (not shown 
here). The experimental data are from Refs. [H. H 13. \2^^ . The data from the PHENIX collaboration denoted by PHENIX 
1 and 2 can be found in Ref. [351 ]. 
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x N^ luon (Q a ) , (2) 



where the first part d dy( ^i pT gives the gluon jet production yield at rapidity y (in pp or AA collisions) computable in 
the kr factorization scheme [13, EH , and the second term N^iuon j g average multiplicity of hadrons in the gluon 
jet with a jet energy Ej et (see Sees. II, III). The symbol <g> indicates a convolution, that is, integrals over variables 
with possible weight factors included. 

The kinematics looks simpler in the center-of-mass of the produced gluon in which two gluons with the mean 
transverse momenta of the order of Q s and the fraction of energy x\ = xi = Q s /\/s collide, producing the gluon 
which moves in the transverse plane with the value of its momentum of the order of Q s . Eq. (|2J up to a possible 
logarithmic correction, can be simply obtained by a dimensionality argument based on the CGC picture in which the 
multiplicity of the gluon jets is proportional to u s Q 2 s where <j s is the effective area of interaction [27]. Notice that the 
typical transverse momentum in Eq. ([TJ is of the order of the saturation scale Q s . 

The crucial ingredient which is essential to explain the different energy dependence of the m ultip licity in pp from AA 
collisions originates from the experimental data for jet production in e + e~ annihilation [29, 36-38], namely N^ luon 
is almost constant at energies of the gluon-jet less than about I GeV but starts to increase with the energy of the 
gluon-jet larger than 1 GeV, see Fig. [2] and Sec. III. For proton-proton collisions, the value of the saturation scale is 
smaller than 1 GeV and consequently N^ luon does not give an additional energy dependence, while for nucleus-nucleus 
scatterings at high energies the saturation scale Q s (A) oc A 1 / 3 Q s (p) (A is the atomic number) is larger than 1 GeV 
and NP luon increases leading to an additional non-negligible power-law energy dependence. This extra contribution 
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FIG. 2: Right: The mean charged hadron multiplicity of unbiased gluon and quark jets in e e~ annihilation, as a 
function of the jet energy . For gluon jet we show experimental data obtained by two different methods, jet boost algorithm and 
subtracting multiplicities in two-jet qq events from three-jet qqg events [36j, 37]. The experimental data for quark jet production 



e annihilation are taken from Ref. [281 ] . The energy behaviour of AT® can be described by E, 



0.6H-0.7 



for En > 0.85 GeV. 



m e 

Left: The ratio of the mean charged particle multiplicities between unbiased gluon and quark jets as a function of scale. Various 
theoretical predictions [3S| based on perturbative QCD (pQCD) are also shown in the plot. The plot in the left panel is taken 
from Ref. 13711 . 



accounts for the gluon-decay effect before hadronization and is missing in the fcr factorization. 

The paper is organized as follows: In Sec. II, we introduce the missing gluon-decay cascade effect in the kt 
factorization. We also show that the observed power-law energy dependence of multiplicity in pp and AA collisions 
at the LHC is fully consistent with the saturation picture by inclusion of the gluon cascade angular-ordering effect. 
In Sec III, we generalize the k t factorization in order to incorporate this effect, and present our numerical results for 
the charged hadron multiplicity both in pp and pA collisions. As a conclusion, in Sec. IV we highlight our main results. 



II. THE ENERGY-DEPENDENCE OF CHARGED HADRON MULTIPLICITY 

The &t factorization [Tt], EH, 113413 includes gluon emissions between the projectile and target, and also gluon 
radiation in the final initial-state from the produced gluons. The &t factorization accounts for the BFKL type gluon 
emissions, namely the parent gluon emits a cascade of gluons with their longitudinal momenta kf being progressively 
smaller while the transverse momenta hn of the parent and emitted gluons are the same. This leads to an angular 
ordering in the cascade shown in Fig. [3] 

p+ > k+ > k+ > ... > fc+, 
p T ~ k T i ~ k T 2--- ~ k Tn , 

ei<e 2 < # 3 < ... < e n . (3) 

However, the other contribution of the gluon decay, in the final initial-state, before hadronization, stems from the 
kinematic region outside the BFKL emission regime where both emitted gluons are collinear to the emitter. In this 
kinematic region, the angle between the gluon (quark) and the decay gluon 0i is small and the main contribution of 
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gluon-decay has an opposite angular ordering to the BFKL type gluon emissions given in Eq. ([3]), see Fig. |3l 

p+ > kf >k+ > ... > k+, 

p T » k T i » k T 2— » k Tn , 

0i>9 2 > e 3 > ... > e n . (4) 

This angular ordering means that a gluon in a fully developed cascade can only emit inside a cone defined by the 
momenta of its first two immediate predecessors, similar to the well-known Chudakov effect in QED (44| . It is well- 
established fact from jet observables (especially at small momentum fractions z = Phadron/ Ej et ) that soft and collinear 
logarithms summed by the Modified Leading Logarithmic A ppr oximation (MLLA), together with angular ordering 
reproduces the most important features of QCD cascade [36l - l38l . |45| . This combined with the Local Parton-Hadron 
duality (LPHD) also gives quantitative predictions for hadron multiplicity and spectra in e + e~ and ep collisions over 
the whole momentum range down to momenta of a few hundred MeV [HI, HU . The MLLA contains systematically 
next-to-leading logarithmic corrections and incorporates single and double-logarithmic effects in the development of 
parton cascades [H, |45[ . 

One should note that although the MLLA angular-ordering kinematic region is quite important at the gluon (quark) 
decay stage, but it does not lead to the large double log contribution in the total cross-section and it contributes to the 
self-energy of the quark and, therefore, to the running QCD coupling in the case that we integrate over all produced 
gluons. Therefore, this kinematic region is not included in the kx factorization formula and has to be considered 
separately, namely, summing of these double log for the gluon jet decay leads to the extra term N^ luon in Eq. (fT]) 
which can be calculated within the MLLA approach. On the other hand, the effect of propagation and interaction of 
the produced jet in the gluonic medium with the BFKL angular-ordering and its saturation effect have been already 
taken into account in the fcx factorization. It is well-known that the gluon decay probability can be factorized from the 
rest of cross-section in e + e~ — > qqg reaction (45[. This is the essence of the factorization given in Eq. (fT]). Therefore, 
one may extract information about the gluon-decay stage in the MLLA region from gluon jet data in e + e~ collisions. 

In order to verify Eq. @, we need to know N^ luon . As we already mentioned, one may calculate the charged 
hadron multiplicities in the gluon jet N^ luon (Ej et ) in pQCD within the MLLA scheme [H, HH- In order to obtain 
the energy dependence of the function ^Giuon^ we uge directly experimental data for N^ luon (Ej et ) in the e + e~ 
annihilation. Unfortunately, such data are limited to high gluon-jet energy Ej et > 5 GeV [H, [37]], see Fig. [3J For 
lower energy Ej et < 5 GeV, we construct the hadron multiplicity of gluon-jet from the corresponding multiplicity of 
quark-jet pjQ uark where we have experimental data [29| . It is well-known that in the double log approximation, the 
ratio of the multiplicities in quark and gluon jets is equal to N^ luon jj^Q uark = Ca/Cf = 9/2 [38|, H^]. However, the 
higher order perturbative corrections significantly suppress this ratio at low energies of the jet making it close to one 
[H IMS, 53, see Fig. H (left panel). One can see from Fig. [2] that N% luon is constant at about E jet < 1 GeV and it 
grows as a power of Ej et at higher energies. From the available e + e~ collisions data shown in Fig. [21 we found that 
the energy-dependence of the mean charged particle multiplicity of gluon-jet can be approximately described by 

(N^ luon ) cx E s jet , with 6 = 0.6-^0.7 for E jet > 0.85 1 GeV. (5) 

It is essential to stress again that such behaviour also follows from the theoretical estimates in the next-to-next-to- 
next-to-leading order (3NLO) pQCD [115153 in the MLLA scheme [M Si- 
It is well-known that the saturation scale has the following energy (or x) behaviour [Id, [ljl I48l - l5l1 | 
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where the saturation scale Qq is fixed at an initial value Xq. We assume that the typical energy of the gluon jet Ej et 
is of the order of average saturation scale. Now, using Eq. ([2]) and Eqs. (15161) we obtain, 

Ql oc fl V2 = s o.H ; (7) 

Ql X (E jet OC Q s )°- 65 OC s V2+0.65xA/4 = s 0.145 5 (g) 

where we assumed that the saturation scale for pp collisions is Q s < 1 GeV and for AA collisions we have Q s > 1 
GeV. In the above and the following we take for the parameter 5, the average value S — 0.65 from Eq. ((S|) and Fig. [5] 
In Eq. (J7]), the average value of A = 0.11 in the the effective saturation scale for pp collisions can be obtained from 
kt factorization results given in Ref. ill] or by a fit to the available data for non-singlet diffractive inclusive hadron 
production in pp collisions shown in Fig. [TJ Then, the power-law behaviour given in Eq. (|5]) for AA collisions comes 
naturally without any extra freedom. In Fig. [TJ we show that the energy power-law scaling given in Eqs. (17181) leads 
to a very good description of experimental data both in pp and AA collisions, including the recent ALICE data in 
A A collisions at 2.76 TeV. 



dNh 
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III. CHARGED HADRON MULTIPLICITY IN THE IMPROVED fc t FACTORIZATION 



In this section, we shall investigate how the saturation model predictions based on the k t factorization [17], [l8j will 
change by the inclusion of the angular ordering effect in the gluon-jet decay cascade. Motivated by previous sections, 
we postulate that the missing effect of gluon-jet decay cascade can be effectively incorporated into the kt factorization 
in the following way, 

(JM, C f J -Gluon 

— ± (AA or pp)= - d 2 PT h[rf\ (AA or pp) M^ luon (Q s ), (9) 

drj a s J dy d z p T 

where h[rj\ is the Jacobin transformation between y and r\ |27| . The impact-parameter dependence of the formulation 
allows to calculate the average area of interaction a s via the geometrical scaling property [IJJ. The gluon jet cross- 
section in AA (or pp) collisions can be obtained from 17] 

da^(y;p T ;B) 2C F a s (p T ) f B2 ^ f rf2 ^ 2 _, ^ V|Ag, (gij r T - b) V^,, (s 2 ; r T ; &_) ^ 



dyd 2 p T (27r) 4 J Bl J ' P%Oi s {QA,p{xi]b))a s (Q A ,p{x2-,b-)) 

where x\ : 2 = (pt / \/s)e ±y , pt and y are the transverse-momentum and rapidity of the produced gluon jet. The 
vector B is the impact parameter between the center of two nuclei (or two hadrons in the case of pp collisions) , b and 
= b — B are the impact parameter between the interacting nucleons with respect to the center of two nuclei (or 
hadrons). A given centrality bin corresponds to a range of the impact-parameter B g [B\,B2] of the collisions. We 
extended the /cy-factorization by introducing a running strong-coupling a s [ll|. In the above, the amplitude N][ p is 
defined as 

NZ p (x % - r T ; b) = 2N A , P (x i; r T ; b) - N 2 AtP (xi;r T ;b) , (11) 

where Na, p (xi]rT',b) is the dipole- nucleus (for index A) or dipole-proton (for index p) forward scattering amplitude 

with tt and b being the transverse dipole size and the impact parameter of the scattering, respectively. 

In Eq. j\f Gluon j g ^ e average hadron multiplicity in the gluon-jet decay in the MLLA region and can be obtained 
from experimental data in e + e~ reactions (29l. l3q. l37j . Using Eq. ([5]) and assuming that typical transverse momentum 
of the gluon-jet is approximately equal to the average saturation scale Q A p at a given centrality and kinematics, we 
define 

fo \ 065 — 
- AW for Qa, p > 0.85 GeV; 

^ luon (Q A , P )=Co{ V J (12) 



with a notation, 



:a,p 



1 for Q A p < 0.85, 



' Q\ p (x u b) + Q\ p (a: 2 ,6-) y /2 
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The normalization factor Co in Eq. (|12p can be absorbed into the parameter C in Eq. §§§ which relates the produced 
gluons to the final-state hadrons based on the Local Patron- Hadron Duality principle |38| , assuming that the final-state 
hadronization is a soft process and cannot change the direction of the emitted gluon-jet further. 

The impact-parameter dependence in the kt factorization is not trivial and a prior is not obvious if it can be 
factorized. Here we are interested to study the effect of new J\f^ luon term in the kt factorization Eq. ©. To this 
end, we employ the b-CGC saturation model [52j which gives a good description of inclusive hadron production in pp 
collisions at the LHC [ll[ . In this model the size of proton naturally changes with energy [ll[ . This model effectively 
incorporates all known saturation properties driven by the small-x non-linear evolution equations [5l| including the 
impact-parameter dependence of the dipole amplitude This model describes very well the HERA DIS data at 
small-x (5ll-l53| and direct-photon production (55|. The extension of this model for the case of nuclear target was 
introduced in Ref. 15] which also gives a good description of RHIC multiplicity data. The dipole- nucleon forward 
scattering amplitude in the b-CGC model |52[ is defined as, 



N p (x; r; b) 



AT (!) 2 ^ + ^M*)) forZ<2; 
1 - exp {-A In 2 {BZ)) for Z > 2, 



(14) 



where we defined Z = r Q p {x; b) , Y = \n{l/x) and K — x" {is) / x' {is) where x 1S the LO BFKL characteristic function. 
The parameters A and B are determined uniquely from the matching of N p and its logarithmic derivatives at Z = 2. 
The proton saturation scale is given by 



) p {x;b') = ( — ) 2 exp 



7/2 



4(1 - lcr )B, 



CGC 



(15) 



Based on the universality of the saturation in the CGC framework, the corresponding dipole-nucleus dipole amplitude 
Na can be obtained from Eq. (fT4|) by only replacing the proton saturation scale by that of the nucleus, 



Qa (x;b) 



d 2 b' T A (b-b') Q 2 p {x;b'), 



(16) 



where Ta {B) denotes the nuclear thickness. The above definition leads to Q\ ~ Q^A 1 ^ 3 which is consistent with 
basic idea of saturation 0, @, [56j]. We use for the nuclear thickness the Wood-Saxon parametrization [57|- The 
parameters A, 7 cr ,iVo, xn and Bqgc a re obtained from a fit to the DIS data at low Bjorken-a; x < 0.01 with a very 
good x 2 /d.o.f. = 0.92 (52|. 

A. Numerical results and discussion 



The number and density of participant at different centralities are calculated based on the Glauber formalism 
assuming cr^ eZ = 64.8,58.5 and 42 mb for yfi = 5.5,2.76 and 0.2 TeV, respectively [H]. Following Refs. [H 
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in order to regularize the divergence of the k t factorization, we introduce a gluon-jet mass m,j e t ■ We are now ready 
to confront the improved k t factorization Eq. ([§]) with experimental data. First, notice that the nuclear saturation 
scale defined via Eq. (fT6|) can be in principle different with exact one upto a factor of the order of one. A change of 
Qa — > 1 -h 1.5 Qa brings about — 7% increase in the hadron multiplicity obtained from Eq. ((9]) at high energies. We 
have only two free parameters, the pre-factor C and the gluon jet-mass rrijet which are determined at low-energy for 
a fixed centrality. The main source of uncertainties in our approach is due to the assumption that gluon-jet mas mj et 
and the normalization pre-factor C do not change with energy, rapidity and centrality. Unfortunately, due to limited 
available data in AA collisions at various high energies we cannot verify if this assumption is correct and therefore we 
should take into account possible uncertainties associated with this assumption. We use RHIC data at y/s = 200 GeV 
around midrapidity to fix our only free parameters rrij e t and C. Unfortunately, as it is obvious from Figs. [1] and [4j the 
experimental errors in the data points taken for fixing these unknown parameters is rather large. We checked that in 
the case of AA collisions, mj et 0.12 -j- 0.14 GeV is consistent with RHIC data within error bars. The experimental 
errors in the charged hadron multiplicity in Au+Au collisions at ^fs — 200 GeV (shown in Fig. [1} may induce an 
uncertainty as large as 7 4- 9% in the value of the parameter C. 

In Fig.Q] we show the energy dependence of the charged hadron multiplicity at midrapidity (labeled with saturation) 
obtained from Eq. © both for pp and AA collisions. The proton saturation scale Q p in the b-CGC model Eq. (fT"5)) 
is rather small and varies very slowly with energy, e. g. for central collisions and midrapidity at y/s = 14 TeV 
and pt = 1 GeV, we have Q p = 0.8 GeV. Therefore, in the case of pp collisions for our interested range of energy 
considered in this paper, the contribution of Nff uon term in the improved k t factorization Eq. ^ is negligible and the 



charge hadron multiplicity obtained via Eq. (0) , shown in Fig. [TJ coincides with the results given in Ref. llj without 
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FIG. 4: Pseudo-rapidity distribution of charged particles produced in Au-Au and Pb-Pb central collisions at RHIC yfs = 130, 200 
GeV and the LHC energies yfs = 2.75, 5.5 TeV. The experimental data are from the PHOBOS [6(| and the ALICE collaboration 



the presence of J\f^ luon term. However, in the case of AA collisions, the nuclear saturation scale defined by Eq. (fTBj) 
can be Qa > 0.85 GeV and consequently J\f^ luon term in the improved kt factorization Eq. ([9]) is important. The 
inclusion of gluon-decay angular-ordering effect via J\f^ luon in Eq. (jSJ does not noticeably affect our prescription at 
RHIC due to our freedom in fitting m je t and C parameters to the same data (at yfs = 200 GeV) while it increases the 
charged hadron multiplicity about 20 — 25% at the LHC energies in AA collisions. Notice that the impact-parameter 
dependence of condition given in Eq. (fT2]) limits the contribution oiN^ luon term at various energies and centralities. 
Overall, the improved k t factorization results Eq. @ shown in Fig. [T] agree very well with both pp and AA data at 
the LHC and also RHIC, including the recent ALICE data for AA collisions at 2.76 TeV. 

In Fig.[4l we show pseudo-rapidity dependence at RHIC energies yfs = 130 and 200 GeV in — 6% Au+Au collisions, 
and also for the LHC energies yfs = 2.76 and 5.5 TeV in — 5% Pb+Pb collisions. In our calculation, the number of 
participant at 5.5 TeV for — 5% centrality is approximately N par = 385. Our prediction for dN^A/dr) obtained from 
Eq. ^ at midrapidity for — 5% Pb+Pb collisions at 5.5 TeV is 1897 ± 133. It is seen in Fig. [3] that as the energy 
increases the peak of rapidity distribution at forward (backward) becomes more pronounced due to the saturation 
effect. This effect has been also observed in Refs. [IU[6l[ m the case of pp collisions. 



In Fig. [5] (right), we show the scaled pseudo-rapidity density (2/N par )(dNAA/ 'drf) at midrapidity where N par is 
the number of participant for a given centrality. The recent ALICE data [|[ at 2.76 TeV A A collisions reveals 
interesting scaling property, namely (2/N par )(dNAA/dr]) at different energies have the same N par dependence upto a 
normalization factor. In Fig. [5] (right), we show that the saturation results obtained via Eq. ([9]) for yfs = 0.2,2.76 
and 5.5 TeV, indeed follow this scaling behaviour. One can observe similar scaling property already at RHIC 62], 
namely dNAA/drj at fixed energy but different centralities falls into a single curve upto a normalization factor, see 
Fig. [5] (left). Both scaling properties shown in Fig. [5] can be easily understood within the CGC picture and follows 
from simple Eq. @. We expect that the centrality-scaling for dNAA/dr] at a fixed energy will be also valid at the 
LHC. Notice that the logarithmic correction due to the running strong-coupling in the kt factorization is not shown in 
the simple Eq. and is important. The curves shown in Fig. [5] are results of full calculation and includes this effect. 
We predict that {2 / N par ){dN aa/ d-q) for 5.5 TeV AA collisions at midrapidity to be about 87±o 06 times bigger than 
the corresponding one at 2.76 TeV, see Fig. [5] (right). 

Finally, one may wonder if a convolution of the fragmentation function in Eq. ([5]), can have the same effect as 
incorporating the gluon-decay cascade effect via _A/"G/uon^ First, one should note that the main contribution of the 
kt factorization for the multiplicity comes from small pr < 1.5 GeV where the fragmentation functions are not 
reliable. On the same line, the fragmentation function is based on a different factorization and QCD evolution 
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FIG. 5: Right: The scaled pseudo-rapidity density as a function of number of participant at midrapidity for AA collisions 
at 0.2,2.76 and 5.5 TeV. Left: The pseudo-rapidity distribution at RHIC 0.2 TeV at different centralities. Both theoretical 
predictions and experimental data show gluon saturation-driven scaling property. We also show in the plots, the corresponding 
normalization product factors. The experimental data are from the PHOBOS UlEl and ALICE [| collaboration. 



equation and its universality is also questionable in the kt factorization approach. Moreover, here we were mostly 
interested to understand the role of initial-state effects in the hadron productions in pp and AA collisions. Therefore, 
we generalized the factorization Eq. ([9]) in order to incorporate the missing initial-state effect due to the gluon-decay 
cascade. We then assumed that the final hadronization is a soft process and will not change the direction of gluon 
decays. This was also motivated by the fact that the MLLA scheme [HI, HU combined with the LPHD principle 
provides a good description of data in e + e~ and ep collisions upto very small pr [H|. 



IV. CONCLUSION 



We showed that the basic energy power-law behaviour given in Eqs. (|7l8p for pp and AA collisions is in accordance 
with the saturation/CGC picture. We showed that the gluon-jet angular-ordering at the decay stage induces extra 
energy-dependence in the case that the saturation scale is large. This contribution has been neglected in previous 
kt factorization based studies. This effect is important for AA collisions where the saturation scale is larger and 
gives rise to an extra contribution about 20 — 25% to the multiplicity in AA collisions at the LHC. This explains the 
observed different energy power-law behaviour of charged hadron multiplicity in AA and pp collisions at the LHC. The 
scaling properties of multiplicity at different energies and centralities shown in Fig. [3] and also the energy power-law 
behaviour of the multiplicity in pp and ^4^4 collisions shown in Fig. [TJ all indicate that the saturation picture and the 
CGC scenario provides a unique and efficient way of describing various experimental data. 

After our paper, there has been recently an interesting paper by Lappi 63] on the same line. It is argued by the 
author of Ref. [(33| that the ratio of (pt) / \J {dN / drj) / St in our approach is in conflict with the observed experimental 
data at RHIC. Here we would like to point out that this is not the case indeed. 

It is constructive first to recall the argument of Ref. [63] . Let us assume that based on the LPHD picture one gluon 
produces n charged hadrons after fragmentation. Then based on only dimensionality argument we have, 

(pr) ~ Qa/n, 
J_dN ^ nQ2 
St dr] 

where (pt) is the average transverse momentum and St is the overlap area between the colliding nuclei in the transverse 
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FIG. 6: Right: The ratio (pt) / \/(dN w /drj)/ St at various centralities for Au+Au collisions at 0.2 TeV. The data was constructed 
from three experimental measurements, the average transverse momentum (pt) of tv~ , dN n /drj and St [6J]. Left: Average 
transverse momenta as a function of y (dN n / drf) / a s for Au+Au collisions at 0.2 TeV. The experimental data are from the 
STAR collaboration [Hj]. 



plane. From above we have (pt) / y 1 (dN/drf) / <r s ~ Therefore in our approach, the ratio (pt) j y (dN/drf)/a s 

decreases for more central collisions in contrast to the KLN type approach [2?], [63|. First, notice that in our approach 
we have n ~ jyGiuon f or > q ^ GeV corresponding to the excess of charged hadron production in the 
presence of jet-decay effects, see Eqs. (|2I12I) and Fig. [2j In Fig. [6] (right) we show the experimental data from the 
STAR collaboration [64J for (pt) / y 1 '(dN* / 'drf) / a s as a function of centrality. It is seen from Fig. [6] that the ratio 
(pt)/\/ (dN / drf) / St at y/s — 0.2 TeV Au+Au collisions decreases for more central collisions. We expect that based 
on our model, this ratio should further suppresses at the LHC in the c entral A A col lisions. 

In the KLN type approaches we have (px) ~ x where x ~ y 1 (dN / drf) / St while in our approach we have 
(pt) ~ x 264 for the case that the saturation scale is Q s > 0.85 -j- 1 GeV. This follows from the fact we have 

x - Q S \J N^ luon (Q s ) and (p T ) - Q s /N^ luon (Q s ) at large saturation scale, see Eqs. (|2I5I12|) . In Fig. [6] (left) we 

show the average transverse momenta as a function of y 1 (dN/drf) / St- The STAR collaboration [64| has found that 
the experimental data for the charged pion at different centralities can be described by (px) ~ Po + 0.07a: where the 
constant po = 0.29 GeV was obtained from a fit and may be interpreted as primordial transverse momentum. It seems 
however that the obtained value of po is rather large (and the coefficient behind x is abnormally small). Such a rather 
large value for po may be in contradiction with the notion of asymptotic deconfinement for a dense system, namely, 
the confinement radius increases with density {1.5, 65]. In Fig. [51 we show that a fit driven by our approach prefers 
a smaller primordial transverse momentum of about pion mass po ~ 0.14 GeV (or even po ~ 0) and it reasonably 
describes the same data within the error bars. Notice that at small multiplicity for very peripheral collisions our fit 
and entire saturation formulation is questionable. 

Finally, we should stress that the jet-decay effects brings rather small extra contribution which requires a more 
careful analysis than only a naive dimensionality argument. We showed in this paper, this extra contribution is 
important when the saturation scale is large and that is in accordance with the existing experimental data in various 
reactions. Our main concern in this paper was the combine description of proton-proton and nuclear data in a contrast 
to the KLN type approaches that deal mostly with nuclear reactions. 
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